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Abstract 

The stability of two quantum counter-streaming electron beams is investigated within the quan- 
tum plasma fluid equations for arbitrarily oriented wave vectors k. The analysis reveals that 
the two quantum two-stream unstable branches are indeed connected by a continuum of unstable 
modes with oblique wave vectors. Using the longitudinal approximation, the stability domain for 
any k is analytically explained, together with the growth rate. 
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I. INTRODUCTION 



The theory of quantum plasma is current 
applications ranging from micro-electronic 



y rapidly expanding result of its potential 
l| to laser-plasma interaction 

as 

or astro- 



physics 

HQ 

. Since the derivation of the quantum fluid equations p, |7| , many well known 
classical plasmas effects including plasmas oscillations with or without magnetic field js, 9|, 



plasma instabilities p, [loll2|. solitons and vortices 13|, or turbulence 14|], have been re- 
examined in their quantum version. 

Beyond their potential applications, many of these problems present a fundamental in- 
terest as they form the basis for the general theory of quantum plasmas. Among them, the 
quantum version of the two-stream instability is endorsing the role its classical counterpart 
plays with classical plasma: that of a fundamental process worthy of investigation in itself, 
and included in every plasma physics textbook. Since the first works on the quantum two- 
stream instability, it has been recognized that one of the most striking difference between 
the quantum and the classical treatment of this process is that the quantum version displays 
two unstable branches instead of only one for the classical case 61. Kinetic treatment are 

n 

now available showing that thermal spread can suppress both quantum branches [15j, as 
they do with the single unstable branch in the classical case. 

While unstable quantum modes in counter-streaming beams systems have been assessed 
for modes with wave vector aligned (two-stream) or perpendicular (filamentation 
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ill) 



with the streams, intermediate orientations have not been considered yet. The intent of 
this paper is precisely to fill this gap. As a result, an interesting conclusion is reached: the 
two branches of the quantum two-stream instability are indeed connected with each other. 
By letting the wave vector acquire a perpendicular component, we found a continuum of 
unstable modes bridging the two quantum branches across the k space. 



II. DISPERSION EQUATION FOR ARBITRARY MODES 

We here derive the dispersion equation for arbitrarily oriented wave vector. For the sake 
of simplicity, we deal with the generic system formed by two cold counter-streaming electron 
beams of equal density no/2 and opposite velocities ±Vqz. Ions form a fixed neutralizing 
background with density rii = n and there's no background magnetic field. The dispersion 



2 



equation is derived from Maxwell's equations and from the quantum fluid equations for the 
two electrons species, 



dt 
and, 



fc + V • (n±v±) = 0, (1) 



+ ( v± • V)v ± = — I E + -± + — j V — ^ , (2) 



<9t m \ c J 2m 2 \ y/n± J 

where — q < and m are the electron charge and mass respectively. The subscripts "+" 
and "-" stand here for the electron population going to the positive and negative directions 
respectively. These equations are now linearized perturbing every quantity with a term 
proportional to exp(ik • r — iut), where k = (k x ,0,k z ). The linearized mass conservation 
equation yields straightforwardly, 

k- v±i f . 

n±i = n± , (3) 

LO — k • V±o 

while the linearized quantum Euler equation reads, 

vi \ 9 ( , v ±0 X BA .h 2 k 2 n ±1 

H k • v ±o - w)v±i = (Ei + )- l -r^2 k - ( 4 ) 

m \ c J Am z n±o 

While flow aligned unstable modes (k x = 0) are exactly electrostatic with Bx = 0, such is 
not the case for a general orientation of the wave vector. This is why the exact calculation 
requires a Lorentz force term. Although some approximations will be used in the sequel to 
analyze the problem, we first outline the exact one. 

The linearized equation ([3H1]) allow to express the perturbed velocities v±i in terms of 
the perturbed fields Ei and Bx- The first order magnetic field is then eliminated through 
the Maxwell- Faraday equation Bx = (c/w)k x Ex so that the first order current Jx = 
~Q S±( n ±o v ±i + ^±i v ±o), can eventually be expressed in terms of Ex only. Combining now 
Maxwell-Faraday and Maxwell-Ampere equations gives, 

(? Aiir 

— kx (kx Ex) + Ex + Jx = 0, (5) 

from which the formal dispersion relation is finally derived, 

■kx(kxE 1 ) + e(k,u)-E 1 = 0, (6) 



c 2 



4«7T 



where e(k, u) = 1 + — Jx is the dielectric tensor. The dispersion equation derived from 
Eq. (jSJ) is very involved and will not be reproduced here. Instead, we now introduce the 








FIG. 1: Growth rate d in terms of K for vq/c = 0.1 and H = 0.6. Top: Exact calculation. Bottom: 
Longitudinal approximation. 



very simple longitudinal approximation which, in the classical case, has been 
over the full wave-vectors spectrum except near the perpendicular direction 
dispersion equation for longitudinal waves reads [181 ]. 

47rg 2 n /2 f /±(v)d 3 v 



bund worthy 
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13| . The 



0=1-5: 



m 



(w-k-v) 2 -(^ 2 /2m) 



(7) 



where f± is the distribution function of each beam. With /±(v) = 5(v x )5(v y )5(v z =F Vq), the 
result reads e(f2) = 1 — F(Q) = 0, where 

1 1 



+ 



(O - K z y - \H* [Kl + Klf (fi + K z f - \H* {Kl + K*) 



in terms of the non-dimensional parameters, 



K 



kw 



n 



H 



hUJ n 



u p Up mv, 



2 ■ 



(8) 



(9) 



where u p = 4irnoq 2 /m is the electron plasma frequency (note that the plasma frequency 
is here defined for a density Uq different from the beam electronic densities n /2). The 




FIG. 2: Growth rate 5 as a function of K for H = 0.6 and 9 = tt/6. 

rescaled growth rate of the instability S in terms of the rescaled wave vector K is plotted on 
Figure [1] without any approximation (top), and according to the longitudinal approximation 
(bottom). As expected, the exact calculation displays unstable filamentation modes for 
K z = which are bypassed by the approximation. Apart for the quasi normal wave vectors, 
the approximation is very good, and even exact for K x = 0. The most striking features here is 
the connection of the two two-stream branches along the flow axis K z . When restricting the 
calculation to this axis, these two branches seem disconnected. After extending the analysis 
for arbitrarily oriented wave vectors, we now see that there is a continuum of unstable modes 
bridging between them. 

III. CONNECTION BETWEEN THE TWO-STREAM MODES AND UNSTABLE 
DOMAIN 

An analytical analysis of the dispersion function for (K X ,K Z ) ^ (0,0) is possible and 
allows for an analytical characterization of the bridge between the two-stream branches. For 
the sake of the discussion, we switch to polar coordinates, 

K x = KsinO, 

K Z = K cos 6, (10) 

and consider 9 = ir/G. Figure [2] pictures the growth rate as a function of K in this direction. 
Figures |3] now display the behavior of the dispersion function for K = 2 (stable), K — 3 
(unstable) and K = 4 (stable). The stability of the system can be explained from the 
respective position of the singularities of the dispersion function. According to Eq. ()8]), 
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these singularities read, 

Q 1 = -Zcos9 --HK 2 , 
2 



Q 2 = -Zcos9 + ^HK 2 , 
Q 3 = Zcos9 -^HK 2 , 

n 4 = Zcose + \hk 2 . (ii) 

For K = 2 in Fig. Ek, we have Q\ < Q 2 < ^3 < ^4- The system is here stable because 
e(Q) > 0, which gives 4 real roots to the 4th degree dispersion equation. For K — 3, the 
ordering of the singularities is still the same, but e(0) < 0, rendering the system unstable. 
Finally, for K = 4, the ordering is Qi < Q 3 < Q 2 < Q4, namely, f2 3 and Q 2 switched their 
positions. The stability in this case is necessary from the behavior of the dispersion function 
near its singularities. When their order is the one of Fig. [3b, i.e. "1,3,2,4", the dispersion 
equation necessarily crosses the Q axis 4 times and the 4th degree dispersion equation has 
therefore only purely real roots. If the ordering is "1,2,3,4", like the one found for K = 2 or 
3, stability relies on the sign of F(0). From Figs. [3b, to [3b, the system is therefore unstable 
until Q 2 = ^3 which gives, 

K* = ^- cos 9, (12) 
H 

defining thus the largest unstable K in a given direction. For 9 = 0, the quantity 2/H is 
exactly the upper limit of the large K instability two-stream branch {fj. For an arbitrary 
angle, the equation above describes a circle of diameter 2/H and center (K x , K z ) = (0, 1/H) 
which fits perfectly as is observed on Figs. [TJ Finally, the expression of K* breaks down 
for 9 ~ 7r/2; according to Eq. ( fl2l . K*(tt/2) = while the fluid quantum theory of the 
filamentation instability gives a finite value for this quantity flO |. 
Turning now the sign of e(0), a little algebra gives e(0) = for, 

V2 



K_{9) = ^—Jcos9 2 - Vcos9*-H 2 , 9 E [0, arccos \fH], 
H 

or, (13) 



V2 



K + {9) = ^—Jcos9 2 + v / cos9 i -H 2 , 9 e [0, arccos Vh). 
H 



From 9 = to arccos VH, K + {9) is monotonically decreasing, and systems with K < K + {9) 
are stable. In the same interval, K_{9) is monotonically increasing, and systems with K < 
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FIG. 3: Dispersion function according to Eq. (|Sj) H = 0.6, 9 = ir/6 and K = 2 (a), K = 3 (b) and 
K=A (c). 
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FIG. 4: Plot of if* (yellow), K- (blue) and (purple) given by Eqs. I[I2HI3|L K- = K + for 
6 = arccos \f~H. Wave vectors between the two frontiers are unstable. 

K_{9) are unstable. The unstable domain defined by Eqs. (TT~2HT~3|) is graphically summarized 
in Figure H] when one recognizes immediately the domain outlined in Figs. [TJ 



IV. GROWTH RATE ANALYSIS 

It turns out that the dispersion equation from Eq. (jSj) is biquadratic and can be solved 
exactly for any parameter set. In terms of Q 2 , the two roots are, 



2 + AKl + H 2 (K 2 X + K 2 ) 2 - 2y/l + AK 2 [2 + H 2 (K 2 + K 2 ) 2 ] 
2 + 4K 2 + H 2 {K 2 + K 2 ) 2 + 2Jl + 4K 2 [2 + H 2 (K* + K*)*\ 



(14) 



Unstable modes clearly arises from Q 2 , since Q 2 , > for any parameters. Of course, the sign 
analysis of Qf yields the unstable domain already explained. 

Analytic results are acessible not only in the electrostatic approximation. Indeed, allowing 
for magnetic field perturbations and following the steps outlined in Section II, the exact 
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growth rate 5 can be found for the filamentation instability (K z = 0), without any further 
restriction, as 



c 2 K 2 x H 2 K^ 2 
1 + — ^ + — - 



v, 







c 2 K 2 



+ Kl[A-H 2 K 2 (l + 



1/2 



l _c 2 Kl_ H 2 K 



On 
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(15) 



From Eq. (1151) it can be shown that sufficiently large wave-numbers such that 



H 2 K 2 



c 2 K 2 

1 + ^4*1 >i 



(16) 

are stable with respect to the filamentation instability. This stabilization effect can exist 
only in the quantum case (H ^ 0). 



V. CONCLUSION 



The essential results of this work are in Fig. 1, where the apparently distinct unstable 
quantum two-stream modes are shown to be connected. This achievement follows after con- 
sidering the two-dimensional k-space and including the Lorentz-force term in the pertinent 
quantum hydrodynamic equations. For arbitrary wave- vector orientation, exact analytic cal- 
culations are possible in the electrostatic approximation, which is fairly reasonable except 
for wave propagation near the axis normal to the beams direction. In this later (filamenta- 
tion) case, the exact growth rate is also acessible. In the general electromagnetic case, the 
quantum filamentation and quantum two-stream instabilities are viewed as a single, unified 
phenomenon. 
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